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Abstract 

We consider the problem of stochastic exit from a planar domain, whose boundary is an 
unstable periodic orbit, and which contains a stable periodic orbit. This problem arises 
when investigating the distribution of noise-induced phase slips between synchronized 
oscillators, or when studying stochastic resonance far from the adiabatic limit. We 
introduce a simple, piecewise linear model equation, for which the distribution of 
first-passage times can be precisely computed. In particular, we obtain a quantitative 
description of the phenomenon of cycling: The distribution of first-passage times 
rotates around the unstable orbit, periodically in the logarithm of the noise intensity, 
and thus does not converge in the zero- noise limit. We compute explicitly the cycling 
profile, which is universal in the sense that in depends only on the product of the 
period of the unstable orbit with its Lyapunov exponent. 
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1 Introduction 

One of the remarkable effects of additive noise on a deterministic dynamical system is 
to enable transitions between otherwise isolated attractors. This phenomenon, known as 
activation in physics and chemistry, is at the origin of the so-called stochastic exit problem. 
Consider a stochastic differential equation (SDE) of the form 

dxt = f{xt)dt + adWt. (1.1) 

Given an attractor A of the deterministic system x = f{x) and a domain T) containing 
A (usually V is taken to be positively invariant under the deterministic flow), the exit 
problem consists in characterizing 

• the distribution of the random time r = inf{t > 0: ^ P} at which paths of the 
SDE, starting in some initial point in T>, leave the domain; 

• the distribution of the exit location x-r on the boundary &D of the domain. 

The best understood situation is the one where / = —VU derives from a potential C/, and 
the attractor A is simply the bottom x* of a potential well. Assume that V contains no 
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other equilibrium point in its interior, and a (non-degenerate) saddle point xi of U on its 
boundary. The expectation of the first-exit time r is then given by 

E{r} = c(a)e[^(-^)-^(-*)l/2-'. (1.2) 

The form of the exponent was already known to Arrhenius [Arr89] , while the value of the 
prefactor c{a), which depends on the curvature of U at x* and xi, was determined by 
Eyring and Kramers in the one-dimensional case [Eyr35, Kra40]. 

Many refinements of this result exist. On the one hand, the theory of large deviations 
allows to compute the exponent in a much more general setting, with a drift term not 
necessarily deriving from a potential [FW98]. In this case, U is replaced by the notion of a 
quasipotential. On the other hand, the computation of the prefactor c{a) in n dimensions 
has been addressed using methods of perturbation theory (see for instance [Aze85, FJ92]). 
Precise results on the relation between the expected first-exit time, capacities and the 
small eigenvalues of the generator of the diffusion have been obtained recently for drift 
coefficients deriving from a potential [BEGK02, BGK02]. 

The distribution of r approaches an exponential one in the small- noise limit [Day83, 
BGK02], and the exit location x-r is strongly concentrated near the saddle xi [FW98]. 
Another related property is that for a fixed time t <^ IE{t}, the probability to leave V 
before time t is exponentially small, of the order t/E{r}. This phenomenon is known as 
metastability, because the state xt may spend extremely long time spans near local minima 
of the potential, while the stationary density has most of its mass concentrated near the 
global potential minima. 

A more difficult situation arises, in two-dimensional systems, when the boundary of 
P is an unstable periodic orbit of the deterministic flow. In that case, large deviations 
theory cannot be applied directly, because dV is a so-called characteristic boundary. The 
situation where dD is an unstable periodic orbit and A an equilibrium point has been 
considered by Day [Day90, Day92] and by Maicr and Stein [MS96]. Day [Day92] proved 
that the distribution of the first-exit location x-r displays a striking behaviour, called 
cycling: As the noise intensity a decreases, the density of Xr rotates around the boundary, 
as a function of |log(T|. Using the concept of the most probable exit path (MPEP) and 
WKB approximations, Maier and Stein also found a periodic dependence on |logcj| for 
the rate of escape per unit time through dT). The intuition behind their estimates is that 
the MPEP spirals geometrically towards &D, until it reaches a neighbourhood of order a 
of dV and escape becomes likely. Where and when this happens depends periodically on 
|logo-|. 

A related situation, which has not yet been analysed in detail, arises when the boundary 
of V is an unstable periodic orbit but ^ is a stable periodic orbit. This situation is 
important in applications: 

1. When studying the dynamics of two coupled, slightly different phase oscillators, the 
onset of synchronization is known to correspond to a saddle-node bifurcation of peri- 
odic orbits. Below a threshold coupling strength, the motion is typically quasiperiodic, 
meaning that the oscillators are not synchronized. Above this threshold, a pair of peri- 
odic orbits of opposite stability appears, the stable one corresponding to a synchronized 
state in which the phase difference between the oscillators may vary periodically but 
is bounded (see for instance [PRKOl]). 

Adding noise to the system will cause phase slips to occur, in which the phase dif- 
ference changes by 27r. Such a phase slip necessarily involves crossing the unstable 
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periodic orbit, so that the determination of tlie distribution of phase shps requires the 
determination of the distribution of first-passage times at an unstable orbit. 
2. The phenomenon of stochastic resonance occurs for instance when a double-well poten- 
tial is forced periodically in time [BSV81]. Noise-induced transitions between potential 
wells are more favourable from the shallower to the deeper well, when the barrier be- 
tween them is lowest. As a result, typical paths of the system will contain a periodic 
component, see for instance [Fox89, WM95, GHJM98]. 

On the mathematically rigorous level, the situation is relatively well understood in the 
adiabatic case, that is for slow forcing, when the paths spend most of the time near 
the bottom of a potential well [FreOO, BG02]. When the forcing is not slow, however, 
it is easy to see that the deterministic system still admits three periodic solutions, 
two stable ones oscillating around the potential wells, and an unstable one oscillating 
around the saddle. Hence the investigation of transition between wells again involves 
the understanding of first-passage times at the unstable orbit. 

One may point out that in the case of synchronization, the phase space has the topology 
of a torus, while for stochastic resonance, it has the topology of a cylinder. However, both 
situations have in common the fact that a path, starting near a stable periodic orbit, has 
to cross an unstable periodic orbit for a transition to become possible. 

In the present work, we focus on the dynamics of the paths up to their first crossing 
of dV. Our aim is to give a precise characterization of the distribution of exit times and 
locations, in particular in the metastable regime. In order to highlight the mechanism 
responsible for cycling, we concentrate here on a simplified, piecewise linear model equa- 
tion, which can be solved exactly to leading order. The general case will be discussed in 
a forthcoming publication [BG03]. 

In our model, we neglect the diffusion in the longitudinal direction, so that the exit 
location is determined by the exit time, modulo the period of the orbit. The main result, 
Theorem 2.3, gives an explicit expression for the density of the first-exit time. If T is the 
period of the unstable orbit, and A its Lyapunov exponent, for a large range of metastable 
times t such that 2|logcr| <^ Xt ^ ^const/a ^ ^j^-g (density is given by 



Here describes the exponential rate of escape provided by large deviations theory. The 
cycling profile P{x) is an explicitly known, universal function, depending only on AT. The 
modcl-dcpcndcnt intrinsic time 9{t), which satisfies 6{t + T) = 6{t) -\- XT, describes the 
"velocity" 1/9' (t) with which the cycling profile rotates around the orbit. 

The regime Xt ^ 2|log(T| is transient, in the sense that paths have not yet reached their 
typical spreading, and thus p+{t) is smaller than (1.3) by a factor e~™"**^ . This 
initial phase is not observed in [MS96], where the authors artificially create a stationary 
regime by reinjecting escaped paths into the attractor. 

Though we do not treat in detail the asymptotic regime t > 6'=°"**/'^", results by 
Day [Day92] imply that the superposition J2k^oP+(^ + ^-^)' which does not take into 
account the winding number k of paths around the unstable orbit, has a similar behaviour 
as (1.3). In fact, we expect (1.3) to hold for times larger than e'^""**/'^ , with an additional 
factor, slowly decaying like exp{— i e"'^^/^'^^}. 




(1.3) 
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2 Model and results 

2.1 Periodic orbits and coordinate systems 

Consider a two-dimensional ordinary differential equation x = f{x) admitting a stable 
periodic orbit enclosed by an unstable one. We assume that the domain S lying between 
the orbits has the topology of an annulus and contains no invariant sets. In that case, 
one can choose polar-like coordinates (r, tp) , such that ip is 27r-periodic and ip > in a 
neighbourhood of S. 

In the deterministic case, it is customary to use as new independent variable, to 
obtain a one-dimensional non-autonomous system 

^=fr{r,<p). (2.1) 

Furthermore, it is possible to choose r in such a way that the stable periodic orbit corre- 
sponds to r = — 1, and the unstable one to r = -|-1. 

Example 2.1. Synchronization 

The onset of synchronization between two weakly coupled phase oscillators can be viewed 
as a saddle-node bifurcation of periodic orbits. A normal- form analysis shows that close 
to the bifurcation point, the dynamics is governed to leading order by an equation of the 
form 

^=r'-ec{^), (2.2) 

where e is the bifurcation parameter. If c{(p) > for all (p, a straightforward analysis 
of the Poincare map shows the existence, for £ > 0, of two periodic orbits r = r+{ip) 
and r = r^{ip) of opposite stability, separated by a distance of order -y/i. The linear 
transformation 



yields the equation 



+ ^^^^^Y-^y (2.3) 

g = i[r+(^)-r_(^)](/-l). (2.4) 

The periodic orbits are now located in y = ±1. Note that the stable (unstable) orbit 
is attracting (repelling) more strongly for those values of (p for which the orbits in the 
original system arc further apart. 

When noise is added to the system, it will in general affect both the transversal r- and 
the longitudinal (/p-direction, so that one is led to analyse the system 

drt = {rt,ipt)dt + ag,. {rt,(pt)dWt, 
dpt = Uin, dt + ag^{ru pt) dWj, 
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possibly witli = 1. Here we consider cr as a small parameter controlling the noise 
intensity, while gr and g^p are fixed functions (of order one). 

Example 2.2. Stochastic resonance 

A classical example of a system showing stochastic resonance consists of an overdamped 
particle in a periodically forced double-well potential, perturbed by additive noise. For a 
Ginzburg-Landau potential, the equation reads 

drt = [n -rl + A cos(27rt/T)] dt + u dWt. (2.6) 

Again, using a Poincare section, one easily proves the existence of one unstable and two 
stable periodic orbits (for any value of the period T). If r-{t) denotes, say, the lower stable 
orbit, and r+(t) the unstable one, a transformation of the form (2.3) yields Equation (2.5) 
with f^p = l,g^p = Q (that is, ip = t), and gr depending only on ip. 

The general case of Equation (2.5) will be discussed in a forthcoming work [BG03]. 
Here we shall concentrate on a simplified model which focuses on the main mechanism 
responsible for the oscillatory behaviour of the first-exit time. 

2.2 Simplified two-level model 

The main motivation to introduce the two-level model is the fact that the distribution 
of first-exit times will be mainly determined by the dynamics near the unstable periodic 
orbit. The dynamics in the remaining phase space can thus be modelled by the simplest 

possible equation, that is, a linear one. More precisely, we will simplify (2.5) by 

• neglecting the term g^, whose effect is a slow diffusion on a time scale l/cr^; 

• neglecting the r-dependence of gr, which is not important near the unstable orbit; 

• replacing fr by a piecewise linear function of r. 

With these approximations, we arrive at the system 

dyt = f{yt,t)dt + ag{t)dWt, (2.7) 

where / and g arc periodic in t, and we may assume that /(±1, t) = for all t. Here time 
t is identified with the angle (p lifted to the real axis. 

In order to avoid certain technical difficulties when dealing with the stochastic process, 
we will actually switch between two linear equations defined in slightly overlapping regions. 
These equations are 

d%- = -a(t)(% +l)dt + ag{t) dWt, 

dy+ = a{t){yt - 1) di + ag{t) dWu ^ ' ^ 

where a{t) and g{t) are T-periodic, positive functions, which are bounded away from zero 
(detailed assumptions will be given in Section 2.3 below). We denote by y^{to,yo) and 
y^(to,yo) the solutions of these equations with initial conditions = yo or = yo, 
respectively. The stable orbit located at y = — 1 and the unstable orbit at y = +1 have 
Lyapunov exponents ^A, where 

A=^^, with a{t)= f a{s)ds. (2.9) 
J- Jo 
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Figure 1 . The process yt is defined by switching between two hnear processes j/j and yf 
each time either the level \ — 5i is crossed from below or the level 1 — ^2 is crossed from 
above. 



The switching between the processes occurs upon reaching levels \ — 5i S (0, 1) from below 
and 1 — ^2 S (0, 1 — 5i) from above (see Figure 1). More precisely, consider the stopping 
times 



= Tyito^yo) = inf{t > to: % (io,yo) > 1 - G [to,oo] 
= T-^ito, yo) = inf{t > to : yti^o^ ^o) < 1 - e [to, oo] 



(2.10) 



Then the process {yt}t^o is defined in the following way: 



yt 



'yr{0,-l) for O^t <Ti = v(0,-l), 
y+(ri, 1 - ,5i) for n ^ t ^ T2 = T^Sn,l - <5i), (2.11) 

,2/r(T2, 1 - for T2 ^ t < T3 = v(t2, 1 - (52), 



and so on. We are interested in determining the distribution of the first-passage time at 
the unstable orbit, namely in the distribution of 

T+ = mf{t>0:yt>+l}. (2.12) 

We use the notation 

p+it) = §^^''-Hr+ ^ t} (2.13) 

for the density of t^, where the superscript in P'^'~^{-} indicates that we start at time 
t = on the stable orbit in y = —1- For brevity, we shall call p+(t) the first-passage 
density of yt to -|-1. 

Note that since is proportional to t (mod T), the first-exit time r+ directly allows 
to determine the first-exit location (pr^. However, r+ actually contains more information 
since it also keeps track of the number of revolutions (or winding number) of the path 
around the stable orbit. 



6 




Figure 2. The periodic functions v*{t), ^^"(t) and 



In the physics hterature, one often considers the rate of escape per unit time, defined 
as minus the time-derivative of the probabihty to be inside the domain D. The difference 
between rate of escape and first-passage density is that the former counts negatively the 
paths which have left D before time t, but returned into P by time t, while these paths 
are not counted by the first-passage density. If the system is symmetric with respect to 
the boundary V, the reflection principles implies that the rate of escape is equal to half 
the first-passage density. 



2.3 Main results 

Four periodic functions will play an important role in the statement of the results. They 
are given by 



1 



W = 1 _ ^-2XT / e-2"(*+^'^) gisY ds, (2.14) 
ft+T 



t+T 



^per 



1 rt+1 

(*) = 3A?n/ e'""+'"^')9(<.)=d«, (2.15) 

where a{t,s) = a{t) — a{s). The first two functions are directly related to the variances 
of y^{to,yo) and y^{to,yo). These are independent of the initial condition yo, and have 
respective values 

a^v^it, to) = [vl^'it) - e-2"(*'*o) vl'^to)] , (2.18) 
a\+{t, to) = [e2"(*'*o) vf'^to) - ^r(*)] • (2-19) 

The function v*{t) allows to determine the qualitative behaviour of v^^{t) and v^^'^{t). 
Indeed, since v^'^{t) satisfies the differential equation v = —2a{t)v + g{t)'^, it is increasing 
whenever it lies below v*{t) and decreasing whenever it lies above. Similarly, v^^^{t) is 
decreasing when it lies below v^{t) and increasing when it lies above (Figure 2). This 
shows in particular that v^'^{t) and v^^^{t) always lie between the maximum v and the 
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minimum v_ of and that they reach their extremal values when crossing v*{t). One 

can also see from the graphical representation that smaller values of the period T lead to 
smaller amplitudes of v^'^it) and v^^'^(t). 
Our main assumptions are the following: 

Hypotheses. 

HI. The functions a{t) and g{t) are twice continuously differentiable, positive, T-periodic 

and bounded away from zero. 
H2. v*{t) has exactly one maximum and one minimum in [0, T), with values v > v_> {). 
H3. There is a constant A > such that v^_^\t) ^ 2v*{t){l-A) and v^'^'(t) ^ 2v*{t){l-A) 

for all t. 
H4. 82/ {2 -82) ^ ^/vJ^. 

H5. has exactly one minimum in [0, T) at a time t = s* , which is quadratic. We 

set 

^ = ^per/^*N inf pP^Vt). (2.20) 

te[o,T) 

Note that Hypothesis H3 amounts to requiring that the most probable exit paths of 
the processes and y'^ can cross at most once the levels 1 — 5i and 1 — 52 (see the 
Appendix). It is always satisfied when T is large enough, because then the variances 
and track v* adiabatically. Is is also satisfied if i; < 2t;(l — A). Hypothesis H4 ensures 
that y^ is dominated by yf near the unstable orbit. We make Hypothesis H5 mainly 
in order to simplify the presentation, it is in fact sufficient to require that the deepest 
minimum of p^'^^{t) be quadratic, which is generically true if v*{t) is nontrivial. 

Under AssTimptions H1-H5, the following result on the first-passage density p+{t) holds 
whenever cr is small enough. 

Theorem 2.3. There exists a ctq > such that for all a ^ gq, the density p+{t) of Tj^ is 
given by 

p+{t) = c{t,a)e-'''/''^\ (2.21) 

where R is defined in (2.20) and the prefactor c{t, a) depends on t and a in the following 

way: 

1. Transient regime: If ^ a{t) < 2|log(T|, then 

1 f Lc""(*h 
c{t,a) ^ const— ^exp< ^ — r> (2.22) 

where L is a positive constant. 

2. Metastable regime: Let Aq = A/(l — A) A {82 — 8i)/8i A 1.-*^ There exists a constant 
l3> such that for t satisfying 2|logcr| ^ a{t) < (t^ e^^o/2<^^ , 

c{t, a) = aCo{s*)e'{t)P (^^^^[^) [l + r{a)] , (2.23) 

where P{x) > is an explicitly known periodic function with period 1, see Equa- 
tions (2.28) and (2.29) below, and 



e{t) = a{t, s*) - I \ogv^'\t) - log 4eF7^, (2.24) 

Z V_^_ [S ) 



^We write o A 6 to denote the minimum of two real numbers a and b. 
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Cq{s*) is a constant given by 



2v*{s*) J ' 



(2.26) 



Si v^a,,(pp-(.*)2) v^j'{s*r/^ 

and the error term satisfies 

r((7) = o(a + ^e-"W). (2.27) 



(7^ 

^2 



3. Asymptotic regime: // a{t) ^ const e^/'^" , most paths will have crossed the un- 
stable orbit as least once, and thus the density decays. 

The periodic function P{x) appearing in (2.23) has the fohowing expressions: 

P(x)= ^ A(Ar(£-x)), with ^(x) = -e-2^exp|--e-2^}, (2.28) 

i=-oo 

which is particularly useful for large T, and the Fourier series 

P{x)= Y: PiQ)e'-''^^, with %) = ^^^rwrfl-^), (2.29) 

q=—oo ^ ' 

where F is the Euler Gamma function. This series converges quickly when T is small. 

Before discussing the implications of this result, let us briefly sketch the proof, the 
details of which are given in Sections 3 to 5. 

A first step is to determine the density V'-l'SjO) of the first-passage time at 1 — 5i, 
when the first switching occurs. We will show in Proposition 4.1 that it is given by 
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V'_(s,0) = -c_(s,0)e-(2-'5i)V2a^^-(s,o)^ ^2.30) 



where a^v-{s,0) is the variance of yj, see (2.18), and the prefactor c_(s,0) does not play 
an important role. Since ?;_(s,0) behaves asymptotically like v^^{s), the first-passage 
times at \ — 5i are sharply concentrated in small neighbourhoods of the local maxima 
of v^^{s). All these maxima correspond to a single point in space, but with a different 
number of revolutions around the stable orbit. In fact, the theory of large deviations 
allows to establish a qualitatively similar behaviour in the general, nonlinear case. 

Next we consider the density t ^ q{t, s) of paths starting at time s in 1 — Si and 
reaching the unstable orbit in +1 at time t. In Sections 4.2 and 4.3, we establish the 
expression 

q{t, s) = -c+{t, s) e-^f/s^'^+C*.^), (2.31) 

where c+(i, s) decays like e"^"*^*''*) and V-^{t, s) = v+{t, s). The exponential decay 

of the prefactor is due to the fact that a large, asymptotically constant fraction of paths 
leave the neighbourhood of the unstable orbit during each period. When computing q{t, s), 
we have to take into account all paths crossing the levels 1 — 5i and 1 — 62 arbitrarily 
often, before reaching +1. This is done with the help of a renewal equation, discussed in 
Section 3. The main contribution, however, comes from paths going directly from 1 — Si 
to +1, without returning to 1 — 62- 



9 



-c{t,a) 




Figure 3. The prefactor c{t, a) of the first-passage density as function of time, for different 
values of a. After a transient regime of duration 2|logcr|/A, the density approaches a 
periodic function of time. Only for times exponentially large in l/cr^ will the density show 
a visible decay. 



For t ^ s, v+{t,s) approaches the function v^'^^{s), which is periodic in the initial 
time s. This means that paths reaching +1 at some fixed time t have left 1 — 5i with 
approximately equal probability near any local maximum oiv^^^{s). This bottleneck ejfect 
is due to the fact that the level 1 — 5i corresponds to a noncharacteristic boundary where 
large-deviation results guarantee concentration of paths, while the unstable orbit causes a 
strong dispersion of paths. 

The first-passage density p+(t) is given by the integral of q{t, s)tp-{s, 0), which can be 
evaluated by the Laplace method, yielding a sum over all minima of the function 

+<^. (2.32) 



v+{t,s) v-{s,0) ' 

For <C s ^ t, this function is close to p^^^{s)'^, and thus has one minimum per period. 
A careful analysis, given in Section 5.1, shows that for t near nT, the feth term of the sum 
is proportional to 



exp{-2(n - k)XT - 2^ e-2("-^)^^} 



2(t2 



(2.33) 



for some 7(t) (the same expression is obtained in [MS96], using WKB approximations). 
This term is the contribution of the paths making k revolutions around the stable orbit 
before reaching 1 — di, and n — k revolutions afterwards. Adding AT to |log(T| will multiply 
0"^ by e"^'^"^, but this only results in a rearrangement of the terms in the sum, without 
changing its value up to small boundary terms. This roughly explains why the prefactor is 
periodic in |log(T|. A more detailed analysis of the sum, given in Section 5.2, is necessary 
to obtain the precise form (2.23). 



2.4 Discussion 
Exponential asymptotics 

The exponential rate occurring in (2.21) is independent of t, meaning that on the level 
of exponential asymptotics, all points on the unstable orbit are reached with the same 
probability. This is in sharp contrast with, say, points on the intermediate level 1 — 6i, 
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Figure 4. The prefactor of the first-passage density as function of llogcrj, for three 
different values oi t = nT. The prefactor is close to a periodic function of llogcrj for 
2|log(T| ^ a{nT) — XnT, and decreases exponentially, with an exponent proportional to 
-const e-"("^) /cr^ for 2|logCT| > a{nT) = XnT. 

which are reached with nonconstant exponential rate (2 — 5i)^/2w_(t, 0). This difference 
is a natural consequence of the fact that the unstable orbit is a characteristic boundary, 
unlike the level 1 — 6i. The rate i?^/2 is the value of the so-called boundary quasipotential 
of the Wentzell-Preidlin theory. 

Note that the expressions (2.21) and (2.23) for the first-passage density are not in 
contradiction with the property 

lim P{e(^'-'^)/2<x2 ^ r+ ^ e^^'+'y^^'} = 1, (2.34) 

holding according to the Wentzell-Freidlin theory for any fixed 5 > 0. This expression 
suggests that first-passage times are concentrated near Kramers' time e'^ ^^"^ , but the fact 
that 6 is independent of a actually allows for a large spreading of the distribution of t+ 
which manifests itself in the behaviour of the prefactor c{t, s). 

Time scales 

Two time scales play a role for the behaviour of the first-passage density p+{t). They 
delimit three time domains, whose boundaries, however, are not particularly sharp. The 
time scale = 2\loga\/X is the time needed for the variances V- and w+ to approach their 
asymptotic values; the effect of the transient phase is still visible, for t ^ t^, in the error 
term 0{e~°'^^^ /<7^)- The metastable time scale tx = Z^*^ , often called Kramers' time, 
measures roughly the time needed for a substantial fraction of paths to reach the unstable 
orbit. There is no sharp transition, however, between metastable and asymptotic regime, 
one rather expects the density p+(t) to decrease geometrically from one period to the next 
one for all times t ^ tr- 

Oscillatory behaviour and cycling 

Consider now the behaviour of the first-passage density p+{t) during one given period 
[nT, (n + l)r] in the metastable regime. The leading term in (2.23) depends periodically 
on t and |log(T|. For fixed noise intensity a, p+{t) is close to a periodic function of t, 
see Figure 3. More surprisingly, for fixed t (which also means a fixed position on the 
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Figure 5. The cycling profile P{x), plotted over three periods, for four different values 
of AT. The normalization is such that the integral of P{x) over one period is equal to 
1/2XT. 

unstable orbit), p+ is proportional to a term depending periodically on |logo"| (Figure 4) 
- this is the phenomenon pointed out in [MS96]. 

It is, however, fundamental to consider the joint dependence of p+ on t and a: The 
leading term in (2.23) can be viewed as a periodic "profile" function of |logcr| — 9{t), 
modulated by the periodic function 0'{t). The maximum of 1 1— > p+{t) moves once around 
the unstable orbit as |logfT| increases by AT: We thus recover the phenomenon of cycling 
described in [Day92, Day96] for the first-exit location, but in addition we obtain here 
the same behaviour for the first-passage time, which keeps track of the winding number 
around the unstable orbit. Also note that the probability of reaching the unstable orbit 
during [riT, {n + 1)T] is proportional to 

,-(,)p(M4^) . ,T P,., d. 4, (2,35) 

where we have set xq = {\loga\ — 9{nT)) / XT . While the peak of p+{t) moves around 
the unstable orbit as a decreases, changing its height periodically, the area below p+{t) 
remains constant. Another way to interpret (2.23) is to consider 0{t) as intrinsic time: 
the first-passage density expressed with respect to the time 9{t) is translated around the 
unstable orbit as a decreases, with constant "velocity", without changing its shape. 

Bottleneck effect 

Paths reach the intermediate level \ — 5i at times concentrated in small windows around 
s* + kT. When approaching the unstable orbit, they are strongly dispersed. As a result, a 
path reaching y = +1 at a fixed time t S [nT, (n -|- 1)T) may have idled along the unstable 
orbit for an almost arbitrarily long time span. The probability that it has come through 
the kth window is proportional to (2.33). When translating this from time to space, all 
windows become superimposed in 93* = (27r/r)s*, meaning that with high probability, all 
paths reaching the unstable orbit have crossed the curve y = 1 — 5i through the same 
small bottleneck. In fact, the same is true for the crossing of any curve bounded away 
from the periodic orbits. We thus recover the well-known fact that the transition between 
the orbits is likely to occur in a small neighbourhood of a fixed trajectory spiraling away 
from the stable orbit, the so-called most probable exit path. But note that the concept of 
most probable exit paths becomes irrelevant for the dynamics close to the unstable orbit. 
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High-frequency limit T ^ 1 

We now examine how the first-passage density changes as a function of the period T. By 
this we mean that the functions a and g are scaled by a factor T, that is, a{t) = aT{t) = 
ai{t/T) and g{t) = gxit) = gi{t/T) for some fixed functions a\ and g\. In particular, 
a{t) = Tai{t/T), so that, for instance, (2.14) becomes 

vPJ'(Ts) = ^ _ e-2rai(s+l,n) ^^(^-,2 (2.36) 

In the limit T — > 0, the functions v^'^ and v^^^ both approach a constant value {gl)/2\, 
where {gf) denotes the average value of gf. The exponential rate approaches the value 
2X[Si + (2 — Si)'^]/{gi) (of course, we cannot actually take the limit T — > 0, because p^^^ 
would also become constant in this limit, so that the discussion is to be understood for 
small but finite T). 

Since |r(l — iy)| decreases exponentially with \y\, P{x) is close to a sinusoid with 
amplitude of order e~™"**/^, so that the cycling profile becomes flat in the high-frequency 

limit, cf. Figure 5. This is related to the fact that many terms, i.e., paths with many 
different winding numbers, contribute to the sum (2.28). In the limit, the cycling velocity 
1/9' (t) behaves like {gl)/Xg{t)'^, so that the first-passage density p+{t) behaves like the 
noise coefficient g{t)'^, independently of a{t). 

Adiabatic limit T » 1 

In the low-frequency limit, v'^^{t) and v^^^{t) both follow adiabatically v*{t), at a distance 
of order 1/T. As a result, the exponential rate approaches [Sf + (2 — 6i)^]/v, which is 
smaller than its value in the limit T — (because 2Xv > {2aiv*) = {gf))- This is due to 
the fact that paths have enough time to probe for the moment when the transition costs 
the least. 

The sum (2.28) is dominated by one term, so that the cycling profile P{x) is sharply 
peaked, cf Figure 5. The first-passage density is thus dominated by paths making a fixed 
number of revolutions around the unstable orbit. The cycling velocity l/0'{t) converges 
to l/a{t). The first-passage density p+(t) thus consists of a peak of height proportional 
to a{t), moving around with velocity l/a{t). 

In spite of the fact that P{x) is sharply peaked, one can show that the first-passage 
density at fixed t is still a monotonously decreasing function of a. This is due to the 
exponentially small factor e~^ /^'^ . 

3 The renewal equation 

In this section, we establish a renewal equation satisfied by the first-passage time r+ at 

-|-1 of the process yt defined in Section 2.2. By restarting the process at the first time 
Ty = T/.(0, — 1) at which the level 1 — (^i is reached, the distribution function of r+ can be 
written as 

P0'-Hr+ ^t} = EO'-i{l{,^^,}P-^'i-^Ur+ ^ t}} 

= f Q(t,s)^_(s,0)ds, (3.1) 
Jo 
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where we have introduced the quantities 

Q{t,s)=¥''^-^'{T+^t} (3.2) 
V^_(5,0) = ^pO.-i{r/^s}. (3.3) 

The first-passage density iIj^(s,0) of to 1 — 5i depends only on the hnear process 
y~; we will discuss its computation in Section 4.1. The function Q{t,s) depends on all 
subsequent switchings of the process yt between and , and therefore we will write it 
as the solution of an integral equation, that will serve as a renewal equation. It is obtained 
by restarting the process each time the level 1 — Ji is reached from below. 

Let denote the stochastic process obtained by killing upon first hitting the level 
+1, and introduce the stopping time 

r*{to) = mf{s >to:y*<l- 62}, (3.4) 

with the convention that r* = 00 if yf either hits level +1 before reaching I — S2 or never 
reaches 1 — 62- 

Proposition 3.1. Q{t,s) satisfies the renewal equation 

Q{t, s) = Pi {t, s) + <3(*> u)K{u, s) du, (3.5) 

where 

P,it,s) = F''^-'^{t+ ^t,T+< T^}, (3.6) 

pu 

K{u,s)= / ^ljAu,v)^ti'"^^)^^^ (3-7) 

J s 

V'/(^,t') = |^F-''-'Hv^«}, (3.8) 
^#^y^s) = -^F-'^'-'^{T*^v}. (3.9) 



Proof: Splitting up the event according to whether > r+ or Tn^ ^ t+, we can write 

Q{t, s) = Pi{t, s) + P^'i-'^iK < ^+ < t}- (3-10) 

Note that Tx^ ^ t+ ^ i if and only if the killed process y^ first reaches 1 — 52 and then 
+1, and both events occur before time t. Therefore, the second term on the right-hand 
side of (3.10) can be written as 

= / / Q{t,u)tpy,{u,v)du tp*{v,s)dv. (3.11) 

Js Jv 

Using Fubini's theorem to interchange the integrals, we obtain the second term in (3.5) . □ 
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Note that since Q{t,t) = 0, the derivatives 

q{t, s) = —Q{t, s), piit, s) = -Pi(t, s) (3.12) 

satisfy the similar renewal equation 

q{t, s) = pi{t, s) + J q(t,u)K(u,s)du. (3.13) 

The renewal equation (3.5) can be solved by iterations. In our case, K will be so small 
that only the first term needs to be computed. The convergence of the iterative method 
can, however, easily be controlled. 

Proposition 3.2. Let Ki{u, s) = K{u, s) and define for n G N 

K{u,v)Kn{v,s)dv. (3.14) 
Then, for any N ^ 1, the solution of the renewal equation (3.5) satisfies 

Q{t,s) = Pi(t,s) + J2 Pi{t,u)Kniu,s)du + RN+i{t,s), (3.15) 

with a remainder satisfying ^ RN{t,s) ^ M^\t — s\^/N\ for some constant M which 
depends neither on N nor ont,s. Hence, Q{t, s) can be written as a series which converges 
uniformly on compact sets. 

Proof: First note that ijj^{t,u) is bounded by some constant M > which implies that 
K{t, s) ^ M for all t, s. By induction, we see that (3.15) holds with 

Rj^{t,s) = Q{t,v)KN{v,s)dv (3.16) 

and that ^ 

^ Kn{v, s) dv ^ M^\t - s\^/N\ (3.17) 

Thus the bounds on RN(t, s) follow from the trivial bounds < Q{t, v) ^ I. □ 

The nth term in the sum in (3.15) is the contribution of those paths which cross the 
levels 1 — (5i and 1 — S2 alternatively n times during the time interval [s,t], before finally 
reaching +1. 



4 The first-passage densities 
4.1 Leaving the stable orbit 

In this section, we ignore all possible switchings between the processes y^ and y^ in 
the definition of yt and focus on the stochastic process y^{to,—l) only. Recall that this 
stochastic process is defined by the SDE 

dyr = ~a{t){yi + \)dt + ag{t) dWt, yi^ = -1, (4.1) 
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t 



Figure 6. The functions -^-(^,0) and uP"(t). 
for some initial time to- Thus it is a Gaussian process, given by 
yt = Vi 



■{to,-l) = -l + a [ e-"(*'^) dT^„ (4.2) 

J to 



where a{t, s) = a{t) — a{s) = a{u) du. At time t, the Gaussian random variable y^ has 
mean —1 and variance a'^V-{t,tQ), where 



v.{t,to)= / e-2"(*'^)<7(5)'ds. (4.3) 



to 

Note that V-{t,to) satisfies the deterministic differential equation 

^t;_(t,to) = -2a{t)v-(t,to)+g{tf, (4.4) 

with initial condition v-{tQ,to) = 0. This equation also admits a periodic solution 

rt+T 



^per 



1 /■*+-' 



Recall that A = a(T)/T is the Lyapunov exponent. These two solutions of (4.4) are 
related by 

v.{t,to) = vnt) - e-2°(*'*o) ^P-(i^)^ (4.6) 

see Figure 6. 

The following result describes the behaviour of the first-passage density ^p-{t,0) of 
y,-(0,-l) to 

Proposition 4.1. Assume that v'^^{t) ^ 2v*{t){l — A) for some constant A > 0. Then 
the first-passage density ip- (t, 0) can be written as 

ij-{t,0) = ic_(t,0)e-''-(*'0)'/2-^, (4.7) 
a 

where 

16 



and 

9itf 



, . 2 - (5i 
c-(t,0) = 



27r 



1 / _ c^—constt^ la^ 



(4.9) 



/or aZ/ 1 < \a{\ - ajl^^) e™"^*^V^' . 

Proof: The process = e"(*)(y^~ + 1) satisfies a stochastic differential equation without 
drift term, 

dz^- = (7 e"(*) g(t) dWt, = 0. (4.10) 

It is a Gaussian process, with mean zero and variance o"^'U_(t,0) = o"^ e^"*-*) t;_(t, 0). The 
first passage of at 1 — (5i corresponds to the first passage of at the time-dependent 
level d-{t) = (2 — 5i)e"(*). Some properties of first-passage densities of such Gaussian 
processes are discussed in the appendix. 

We will apply Corollary A. 5 from the appendix, after checking that ^-(t, 0) and d-{t) 
satisfy the required conditions. (A. 11) can easily be verified by rewriting the condition with 
the help of v^^^ (i^) and recalling the assumption on the relation between t;^'"^(t) and v*{^). 
Note that A in this assumption and A in (A. 11) differ by a constant. Checking (A. 13) 
is even easier, just keep in mind that 7;_(t,0) approaches (i?) and thus is bounded. 
Next, (A. 12) can be established by showing first that there exists a 5o > such that |c(t, s)\ 
as defined in the appendix is bounded above by constat — s as long as t — s ^ Sq. For 
t — s > So, one can show that \c{t, s)\ is at most of order one. Finally, Assumption (A. 22) is 
seen to be satisfied for large enough M3 by comparing again v-{t, 0) with v^^(t). Applying 
the corollary immediately yields (4.7)-(4.9). □ 

4.2 Reaching the unstable orbit 

We now turn to the process 7/^(s, 1 — 5i), defined by the SDK 

dy+ = a{t){y+ -l)dt + ag{t) dWt, y+ = 1 - ^i, (4.11) 
for some initial time s. It is given by 

y+ = y+{s, 1- Si) = 1-6, e"(*'^) +a^*e"(*'") g{u) dWu. (4.12) 

At time t, y^ is Gaussian with variance a^v-\-{t, s), where 

v+{t, s) = e^"^*'") g{uf du. (4.13) 

As y^" is spreading fast as t increases, it is helpful to consider the stochastic process 

z+ = 5i+e-"(*'^)[2/+-l], (4.14) 
which is the solution of the SDE 

dz+ = a e-°^(*'^) g{t) dWt, zf = 0. (4.15) 
It is also Gaussian and has variance a'^v^{t, s), where 

v+{t, s) = e-2«(*.^) v_^_[t, s) = e-2°^("'") g{uf du. (4.16) 
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Note that s i-^ v^{t, s) satisfies the deterministic differential equation 

d 

^v+{t, s) = 2a{s)v+{t, s) - g{sf, (4.17) 
with v+{t,t) = 0. This equation also admits a periodic solution 

vl^is) = £^^e2-(^+^'«) giuf du. (4.18) 

The function v+{t, s) can then be determined by the relation 

v+{t, s) = vl^'is) - e-2°(*'^) vl^'it), (4.19) 

see Figure 6. By the reflection principle, the distribution function of the first-passage 
time f+ of yf' at +1 is given by 

ps,i-5i|~_^ ^ = 2¥''^-^^{y+ ^ 1} = 2$(-^±^), (4.20) 

where ^>(x) = (27r)^^/^ e~"^/^ du denotes the distribution function of the standard 
normal law, and 

i2 2a{t,s) r2 

v+{t,s) v+{t,s) 

The density of f+ can thus be written as 

dt ^^+^*^"aV2^^+(t,.)3/2*' 

Let us now show that the function 



e-p+{t,s)V2a^ . (4.22) 



Pi{t,s) = ^P^'1-^Ht+ ^t,r+< r^} = ^P^.i-'5i{f+ ^t,f+< r^} (4.23) 

behaves in the same way. Recall that t-^ = t-^{s,1 — di) denotes the first-passage time 
of y^{s,l — Si) at the lower level 1 — 52. It is obvious that pi{t,s) is bounded above 
by (4.22), but we want to show that it actually has the same exponential asymptotics, 
and, moreover, almost the same prefactor. In other words, it is not only unlikely that a 
path reaches the unstable orbit, but also the conditional probability that a path returns 
to 1 — (^2 before reaching -|-1, given that it actually reaches the unstable orbit at -|-1, is 
small. Hypothesis H3 plays a crucial role here. 
In a first step, we study the density of r^^. 

Lemma 4.2. The density tp^{u,s) = ■^^^'^~^^{t-\. ^ u} of = t^{s, 1 — Si) satisfies 

^^(n, s) = ic^(n, s) e-P^i.u,s)V2c^^ (4_24) 
where the prefactor c^{u,s) is bounded by a constant times u+(u, s)"^/^ e""^"''^-' and 

v+{u,s) 
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Proof: Recall the definition of the stochastic process = (5i+e-"(«'")[y+-l] from (4.14). 
The first passage of y+ at 1 — 82 coincides with the first passage of z'^ at the time-dependent 
level d^{u) = 5i — ^2 e""^"''') . Therefore, the result follows from Lemma A.l of the ap- 
pendix which we apply with d{u) = d+{u) and v{u) = ?5+(n, s) for fixed s. □ 

Proposition 4.3. Assume that v^^^{u) ^ 2v*{u){l — A) for all u. Then 

pi{t,s) = -c+(t,s)e-''+(*'^)'/2-', (4.26) 
where p+{t,s) is the function defined in (4.21) and 

c+{t,s) = ,,,, - io(e-«(*.^)e— *^o/-^), (4.27) 

V27r v+{t,sYi^ G 

where Aq = A/(1 — A) A (52 — (5i)/(5i A 1. // the condition on is not satisfied, (4.26) 
remains valid, hut the pref actor can he smaller than (4.27) . 

Proof: The main idea of the proof is that the condition on excludes that the most 
probable path going from (s, 1 — 61) to (t, 1) crosses the lower level 1 — 52- The only 
difficulty resides in exploiting this fact in order to obtain the exponentially small error 
bound in (4.27). 

By definition of Pi(i, s), 

Pi(t, s) = P^'i-^i{f+ ^ t} - P^'i-'^i{rx, < f+ ^ t}. (4.28) 

Recall that denotes the first-passage time at level 1 — 82 for the process yY , which is 
obtained by killing upon reaching level -fl. The second term on the right-hand side 
can be written as 



= I ¥''^^-^^{f+^t}^t^{u,s)du. (4.29) 

J s 

Since P*'i-'^2{f+ ^ i} = and V'*(n, s) ^ Vn.('"> «), we have 

^p.,i-5i|^^ <f^^t}^J^ ^F^'^-^^{f+ ^ t}^^{u, s) du. (4.30) 

The first term in the integrand is similar to (4.22), the only difference lying in the initial 
condition. We can write it in the form 

^p«,i-52|^^ ^ = -c+(t,u)e-'^+(*'")'/2a2^ (4 3;^^ 

Ot (J 

where p+{t,uf = 8^/v+{t,u) and c+{t,u) = 525'(t)^ e^^^^*'") /^/2^^v+{t,uf/'^ . The deriva- 
tive pi{t, s) thus satisfies 

cr Lv27r u+Ct, s)-^/^ a 

(4.32) 
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where 

Assume that t — s > 2. (For t — s ^ 2, it suffices to apply the first one of the arguments 

below.) We spHt the integral in (4.32) at ni = s + 1 and U2 = t — 1. The integrals 
over [s,ui] and [^2,^] are easily seen to be exponentially small because x{^) diverges as 
u — > s and u — ^ t (in fact, for small enough a, the integrand is maximal at Ui or U2, 
respectively). In order to bound the remaining integral over [ui,'U2], we only need to find 
a positive lower bound for the function x(^)) valid whenever s ^ u ^ L Using the fact 
that v+{t,u) = e^"('''*)[u+(i, s) —v+{u,s)], it is straightforward to show that 

r{u){l — r{u)) 

where ^ 

,^(u) = ^ e-"("'^) and r(u) = G [0, 11. (4.35) 

£>! v+[t,s) 

We may assume the existence of a constant c* > 1 such that u+(t, u)/v+(t, s) ^ c* for all 
u G [s,t]. The function r{u) being increasing with range [0, 1], we can define u* G {s,t) by 

1 _ r{u*) = e-2"(»*.^) Ml}^ = 1. (4.36) 
Consider first the case u* . Then 



^{^^] i^-r{n)), (4.37) 



and thus 



l — r(u) \Si J v-\-{t,u) \Si J c* \6i 

Ku)-(I-.W)f ^^/^_ 

1 — r{u) c* \di J c*df 

We turn now to the case u ^ u*, in which it suffices to find a lower bound for 9{u) = 
k{u) + r{u) — 1, since x(^) ^ 4:9{u)^ for all u. First note that 9{s) = {S2 — Si) /Si and, 
since (4.36) implies e~"("*'*) ^ 1/c*, we also know that 9{u*) ^ {82 — 5i)/c*5i. Thus if 
9{u) reaches its minimum on the boundary of [s,it*], we are done. Otherwise, the fact 
that 

-2a{u,s) 

9'(u) = ^ k(u)=2v*(u)— (4.39) 

v+{t,s) 

shows that if 9{u) reaches its minimum in {s,t), then 

2v*{u 



g-2a(u,s) 2 



v+{t, u) 



- 1 



where the last estimate holds due to our assumption on v^'^'^{u). Note that if v^'^^{u) is 
not bounded away from 2v*{u), then the minimum of 9 is still positive, but may be of the 
order e~^"(*''^) , which can become very small. □ 



20 



4.3 The renewal kernel 

The kernel K{u, s) involves the functions ipy{u, v) and ip^{v, s) defined in (3.8) and (3.9). 
Wc already know from Lemma 4.2 that the rate associated with ip\^{v, s), which provides 
an upper bound for ip^{v, s), is given by 

PAv,sf = J ^ . 4.41 

Note that p~^{v,s) vanishes at a time v ^ s such that e""'-^'*) = 5i/52, which shows that 
most paths starting at level 1 — (5i at time s will reach the lower level 1 — 62 close to that 
time. 

The same argument as in Proposition 4.1, applied for a different initial condition, 
shows that ip y{u, v) can be written as a~^c y a rate 

, ,2 [(2-(5i)-(2-^2)e-"("'-y 
Py{u,v) =^ r (4.42) 

V- [U, V) 

which is bounded below by {S2 — 5i)^/v. It thus follows from the definition (3.7) of K that 

K{u,s) (con.si/(7)c-(^2-^i)'/2^^'. 
By differentiating (3.15), we find 

00 „t 

qit,s)=pi{t,s) + Y^ / pi{t,u)Kn{u,s)du. (4.43) 

(Note that the sum is converging uniformly on compact sets.) The smallness of K is not 
yet sufficient to ensure the smallness of the sum, relatively to pi{t, s). The following result 
provides a sufficient bound. 

Proposition 4.4. Let p^{t, s) := ^{c+{t, s) \/ 1) e-P+^^'^^^Z"^"^ ^ Pi{t,s)? Assume that 
v+^^{u) ^ 2v*{u){l — A) for all u. Then, the relation 



Pi{t, u)Kn{u,s)du ^ 



const (i-s)^e-™"^*^o^' 



(72 



Piit,s) (4.44) 



holds for all n^l, with Aq = A/(l - A) A {62 - Si)/Si A 1. 

Proof: The main idea consists in comparing ijjy{u,v) to the density il)'y.{u,v) of the 
first-passage time of yf{y^l — 82) at 1 — 5i (recall that tpy(u,v) relates to the process 
y^{v, 1 — (^2))- The exponential rate p'^{u,v) associated with ^p~^{u,v) is given by 

We claim that 

pAu,vf ^ p+{u,vf (4.46) 

for all u V, as a consequence of Hypothesis H4. Denote e""*^"'") by and v/v by k^. On 
the one hand, it follows from the definitions (4.3) and (4.16) of V- and v+ that 

v-{u,v) ^ g{wf dw ^ -^v+{u,v). (4.47) 
Jv s 



^We write o V 6 to denote the maximum of two real numbers a and b. 
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On the other hand, since g{w)'^ = 2a{w)v* (w) and v*{w) G \v,v], we can write 

lv-{u, v) ^ r 2a{w) e-2°(«'"') du; = 1 - = T 2a{w) e-2"("''") dw ^ v). 

^ J V J V H 

(4.48) 

We can thus conclude that 

^^^^^^ 

Using the monotonicity of the ratio in the preceding estimate, we see that Relation (4.46) 
is satisfied provided [(2 — 6i) — (2 — 62)^]^ ^ S2 — for all ^ ^ n, which easily follows 
from Hypothesis H4. 

In order to prove (4.44) for n = 1, we first consider the rate /5^{t,v) associated with 
^v,i-S2_!^f^ ^ compare (4.31). Proceeding as in Proposition 4.3, it is straightforward 
to show that 

p+{t,vf ^ p+{t,uf + p+{u,vf 

^p+{t,uf + p^{u,vf (4.50) 

for all u G [v,t]. This is actually a consequence of the Markov property and the fact that 
yf{v, 1 — 62) has to cross the level 1 — 61 before reaching 1+. Equation (4.50) implies that 

J Pi{t,u)K{u,s)du 

J S J V ^ 

^ const J e-^+(*'^)'/2'^' il;-^{v, s) dv 

^ const ^ e-^°"^*^o/<^' g-p+(t,s)V2a2 ^ ^4 5^^ 

where the last line is obtained by the same argument as in Proposition 4.3, compare (4.32). 
This proves (4.44) for n = 1, and for larger n the result follows easily by induction. □ 

As a direct consequence of Proposition 4.4, whenever t — s < (j2 g const A^/o-^^ 
bound the sum in (4.43) by a geometric series, and conclude that 

-c+{t, s) e-^+(*.-)V2a2 ^ q^t^ 5) ^ ic+(t, s) e-^+(*'^)V2^^^ (4.52) 
a a 

where c+{t, s) = c+{t, s) + 0{{t - s)a-^{l + c+{t, s)) g-c""^*^?/'^') . 

5 Properties of the exit law 

The first-passage law of the process yt at +1 can be expressed as a function of the first- 
passage density ip-{s,0) to 1 — (5i and q{t,s) via (3.1) as 

p+{t) = ^ *> = /* ^)^-(^' 0) (5-1) 

From now on, we assume that t ^ for some sufficiently small constant /? > 0, 

so that the error terms in (4.9) and (4.27) are at most of order h{a) = ae~f^^ol'^" ^ and 
the one in (4.52) is at most of order h{a){l + c+(i, s)). 
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5.1 Estimating the integral by a sum 

In virtue of Proposition 4.1, Proposition 4.3, and (4.52), the integral in (5.1) can be written 



as 



p+(t) = J- Ac+(t,s) + C»(/i((7)[l + c+(t,s)])]c_(s,0)e-[''+(*'^)'+''-(^'°)'l/2'^'ds. (5.2) 
<7 Jo 



The exponent in (5.2) is of the form 

51 ^{2-5,f 



v+(t,s) v-{s,0) 

6j (2-<5-^2 



V 



per 



^ ..per/ N -.Per/r.\^ W-^V 



(s) _ e-2a{t,s) ^Per(^^ ^per^^^ _ ^_2a(s) ^Pe=-(o) 



which we write for ^ s ^ t as 

p(o)(t, sf = pP'^'isf + 0(e-2°^(*'")) + 0(e-2«W) (5.4) 

with a periodic part 

pPer^g^2 _ _]_ _| per, . , (5-5) 

v!^ (s) w_ (s) 

cf. Figure 7. In Hypothesis H5, we assumed, for simpUcity, that pP*^''(s) has a unique 
minimum in the interval [0, T) , at some s* G [0, T) satisfying 

pP^Ys-^) = inf pP"(s). (5.6) 
se[o,T) 

When the time interval [0, t] includes many periods, p^^^ {t, s) will have minima of compa- 
rable depths near all times <^ s* + kT <C t, sec Figure 7. Other minima, which may exist 
for s near and t are much shallower, and thus contribute less to the integral (5.2). The 
integral will be of the order (i/r)c7 e"''''"^**^ /^^^ . On the other hand, when t is not large 
enough, even the deepest minimum of p^^\t, s)^ will be substantially larger than pP^'"(s*)^, 
leading to a value of p+{t) which is orders of magnitude smaller than (t/T)a e"^^""^^*^ ^^^^ . 
The system is then still in its initial transient regime. 

The transition between these regimes occurs when a{t) ~ 2|loga"|. We first show that 
for times such that a{t) is smaller than 2|logcr|, the density p+{t) is much smaller than 
its "asymptotic" value {t/aT)e-P'""'^'*^^/^''\ 

Proposition 5.1. Assume that a{t) ^ 2i^|log(T| for some u < 1. Then 

p+{t) ^ const^e-^'^'''-'^' e-pP-(.*)V2.^, (5.7) 

where L > is a constant independent of a, t and u. 
Proof: From (5.3) we obtain 

pW(i, sf ^ p^-{sf + 5? +(2 - ^1)^5^ e-M^) . (5.8) 

[sy v_ [sy 
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Figure 7. The functions pP™(s)^ and p^'^\ti,sY for four different final times ti < t2 < 



The periodic functions v'^'^[s) and v^"^ [s) being bounded above by f , we have 

where /3i = bivf^ty^/v and /^a = (2 - (5i)t;P''(0)i/2/7;. The right-hand side of (5.9) 
reaches its minimum when s satisfies [32 e~^"^*^ = f3i e~"^*^ , and has value 

pv-^[s* f + 2/3i/32 e-2"W ^ /-(s*)^ + 2fii(32CJ^r (5.10) 

We can now estimate the integral (5.2) by splitting it at times ti = \ and t2 = t — \. (If t ^ 
2, the argument is even simpler.) The integral over [ti, can be bounded by using (5.10), 
while the integral on [0,ti] is small because for a small enough, c_(s, 0) e~''~^*''^^ /^'^ 
reaches its maximum at ti and the remaining factor in the integrand is bounded. On \t2,t\ 
the situation is similar. □ 

We now turn to the case where t > 2|logcr|/A. Assume that t £ [nT, (n + 1)T). If we 
can show that p^^\t,s) has exactly one minimum in [kT, {k + 1)7") for 1 k n, 
and that this minimum is quadratic, the Laplace method will allow us to approximate the 
integral (5.2) by 



-(T.\ L ( ml .,. c+(t,.,)c„(gfc,0) + error term) ^-p''Ht,s,? . (5.11) 



Making this argument precise, we obtain the following result. 

Proposition 5.2. Assume that t G [nT, (n + 1)T) with nXT ^ 2|logcr| and n ^ 4. Then 

P+{t) = lc(.^)4S^5(n,a,t)[l + 0(a)] e-''^^^^^*)^^'^^ (5.12) 
a v\_ [t) '- 

where C{s*) is a constant given by 



Si v^a..(/,p-(.^)2) vi'^'is^r/^ 



2v*{s*) 



(5.13) 
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and S{n, a, t) is the sum 

S{n,a,t) = 2|l^exp|-2(n-A.)AT- 2^[7oe-2'^^^+7(i)e-2(--^)^^]|, (5.14) 



k=l 

with 



-fo = {2-Sife 



-2a(s*) 



v_ 



^per / ,\ 

7(0 = <^?e--(*-^+'^^)^S^. (5.15) 

Proof: We split the integral (5.2) at times ti = kiT and t2 = (n — ki)T, where ki < n/2 
will be chosen in such a way that for ti ^ s ^ t2, p^^\t, s)"^ has a minimum close to s* + kT 
on each interval Ik = [kT, [k + 1)T), while the contributions of the integrals over [0,ti] 
and [t2,t] are negligible. 

Take ki of the form ki = \iy\log a\/2XT] V 2, with a parameter v yet to be chosen. For 
ii ^ s ^ t2, one has e-2°{«) ^ a'' and e-2"(t.s) ^ ^-i^ ^e first show that p^^^t,s)'^ has a 
quadratic minimum in each 1^, ki ^ k < n — ki. For this purpose, we write 

JO)/, „^2^ ^! I ~ '^i)^ /5i6^ 



per 



{s) - eMs,kT) ^Per^g) _g-2a(s,feT)^2' 



where 71 = g-^C^-'^)^^ e-^"^*'"'^) t;^*"' (t) and 72 = e^^^^^^ ?;P''''(0) are at most of order a''. 
Were 71 = 72 = 0, (5.16) would reduce to p^^'^{s)^, which has a unique minimum in /fc 
at s* + kT, the latter being quadratic. Hence, the implicit-function theorem applied to 
dsP^'''\t, s)"^ shows that for sufficiently small 71 and 72, p^^\t,s)'^ has a unique minimum 
in 4 at a time Sk = s* + kT + Oia"). 

Expanding (5.16) into powers of 71 and 72 shows that 

pW(t, skf = p^'^'is*)' + 7(t) e^2(n-fe)AT ^-2kXT +o{a^'), (5.17) 

for the coefficients 70 and ^{t) given in (5.15). Evaluating the integral (5.2) restricted to 
the interval [^1,^2] by the Laplace method yields 




4-7r 

[c+{t, Sk) + 0{h{a)[l + c+{t, Sk)])]c.isk, 0) [1 + Oia^)] 



5..(pP-(..)2)L^--'"- ' -^-^^^^ ■ , ,j ^^^^^^ 



k=ki 



where 5,,pP^''(sfe) = 5,,/)P'=^(s'^) + Oia'') and 

(^1 r7('/'l2 p-2a(t,s*+nT) „-2(n-fc)AT 



c_(sjfc,0) = 



2 -Si 



1 



vr{s*) 2v*{s*) 



*\2 



vP_^'{s*y/^ 



[1 + Oia" + h{a))] . (5.20) 



The error term 0{a^'^) in the exponent (5.17) yields an error term 0{(j'^^'^^^^) in the 
sum (5.18). For this reason, we choose = 3/2, so that all error terms are of order a at 
most. This choice of u is always possible, since 2ki ^ 3|logc7|/(2AT) + 2 is smaller than 



25 



n by our condition on |log(T|. Note that the additive error term 0{h{a)) in (5.19) can 
be incorporated into the sum S(n,a,t), as nh{a) ^ 0(th(a)) ^ 0{a^) ^ 0(a'^)S{n,a,t), 
where the last inequahty will be proved in the next section. 

We now turn to computing a bound for the integral (5.2) restricted to the interval [0, ti]. 
We first consider the case |logcj| ^ 5AT, in which fci ^ 4 and c^^"^*^) = e^^'^'^'*'-^ ^ a^^^. 
With the help of (5.9), the integral can be estimated by exp{— pP'^^(s*)^/2c7^} times its 
maximal value which itself is bounded by 

^exp|-^e-"(*4 =of--i-^exp|--4lV (5-21) 

Hence the first part of the integral is exponentially small compared to the integral itself. 
A similar estimate shows that the integral over [^27^] is small. 

In the case |logo"| < 5AT, the location of the minima has to be estimated with 
more care. In this case, e"^'^"^ ^ cr^/^, so that the same argument as before shows 
that p^^^it^s)^ has a minimum in each of the intervals /i, I1 and /a, and possibly also 
in /q- An examination of dsP^^^ {t, s)'^ shows that these minima are actually located in 
Sk = s* + kT + 0(e-2°(s*+fer))^ ^ ^ 0,1,2,3. We set n = a{s*)/XT, ^ k < 1, and 
distinguish between two cases: 

• If A{k + k)\T ^ 3|logc7|, then p^^\t,Skf = p^^'is^f + 70 e'^'^^^ +0(^3), and the 
integral between and Sk is exponentially small. 

• If 4{k + k)\T < 3|logo"|, then e'~^"(*'=) ^ o"^/^, and the error term only leads to an 
error of order a in the prefactor. 

The intervals Ik, n — A ^ k ^ n — 1, are treated in a similar way. 

The same type of arguments shows that the error made by extending the sum in (5.18) 
to all A; between 1 and n is also negligible. □ 



5.2 Properties of the sum 

To complete the proof of Theorem 2.3, we will show that the sum S{n,a,t), defined 
in (5.14), is close to a periodic function of |log(T|. In this section, we always assume that 
t € [nT, (n + 1)T) and that nXT ^ 2|log a\ as well asn^ 4. To highhght the cr-dependence 
of S, it is convenient to set a = and S{n,a,t) = a^S{n,r),t). Further introducing the 
notations 

^o = -^log7o, m = ~^og-f{t), (5.22) 

and changing the summation index from k to £ = n—k, allows to write the sum in compact 
form as 

n-l 

S{n,ri,t) = ^A{e\T - r] + e{t))B{{n - e)XT - r] + 9o), (5.23) 

with 

A(x) = iexp|-2x-^e-2-|, (5.24) 

B(x) =exp|-^e-2^|. (5.25) 

The function A(x) decays like c^^^ as x ^ +00, and like exp{ — ^ e^'^'l} as x ^ —00. 
It reaches its maximal value c^^ when e~^^ = 2. The function B{x) is monotonously 
increasing. It decays like exp{ — ^ e^'^'} as a; — ^ —00 and approaches 1 as a; — ^ +00. 
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Since 7(t) involves the location of t € [nT, (n + 1)T) relatively to s* + nT (sec (5.15)), 
9{t) is a (right-continuous) periodic saw-tooth function, making a jump of —AT at each 
integer multiple of T. We can thus write it as 



e{t) = e{t) - XT 



(5.26) 



where \_x\ denotes the largest integer smaller than or equal to x, and 9{t) is a continuous 
function given by 

e{t) = a{t, s*) - I log??f ■■(*) - log 4^7^- (5.27) 

Note that e{t + T) = e{t) + AT. 

It is easy to show that S{n, rj + XT, t) — S{n — 2,r],t — 2T) is exponentially small (just 
shift the index of summation and show that the boundary terms are exponentially small). 
The following result gives a more precise characterization of S{n, rj, t) by showing that it 
is actually close to a periodic function of r] — 9{t). 

Proposition 5.3. Assume that nXT ^ 2ri. Then 

Sin,r],t) = S{q,t)[l + 0{a^^)], (5.28) 
where fi = fi{n, rj) = {nXT — 2r])/r] and 

OD 

S{r],t)= A{£XT -r] + e{t)). (5.29) 

e=-oo 

Proof: We split the sum S{n,r],t) into two parts. For ^ £ ^ ra/2, we use the fact that 
1 ^ B{{n-i)XT-ri + eo) ^ 5(|Ar-r? + ^o) = expj-^ e^'?"'*^^} = l-Oia"). (5.30) 

Hence replacing i? by 1 in the sum for these values of i only yields a multiplicative error 
of order 1 — 0{a'^). For £ > n/2, it is obvious that 

A{eXT -n + m) ^ ^ e-2^^^+2'' . (5.31) 
Bounding i? by 1 allows to bound the sum over £ = [n/2] , . . . , n by the geometric series 

n , s n 

A{£XT-r, + m)^^ E = O(a^). (5.32) 

e=\n/2'] e=ln/2] 

Thus the main contribution to S stems from £ G {1, . . . , [n/2j}, and so does the main 
contribution to Yl'^o ^{^^'^ ~ + ^(*))- It remains to check that the contribution of 
negative £ to S" is small. Comparing that sum with an integral shows that it is in fact of 
the order e"'''^*)/^'^ . Finally, replacing 9{t) by 9{t) only results in a shift of the summation 
index. □ 

The function S{r],t) is clearly periodic m rj — 9{t) with period AT. Let us thus write 
S{r],t) = P{{r] — 9{t))/XT), where P{x) > is periodic with period 1. It remains to 
compute the Fourier series of P{x). 
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Proposition 5.4. The periodic function P{x) admits the Fourier series 

oo 

P{x)= ^ P(g)e2-'?-, (5.33) 

g=— oo 

where the qth Fourier coefficient is given in terms of the Euler Gamma function by 
Proof: We have 

P{q) = f P(x)e-2^'«^dx 
^0 

oo „i 

= V / A((^-a;)Ar)e-2'^'«^da; 

£=-oo"^0 

/oo 
A{XTx)e^'''"'''dx. (5.35) 
-oo 

Replacing A by its definition and using the change of variable z = e~^^^^/2 yields the 
result. □ 

Appendix 

Let v{t) be continuously difFerentiable on [0, oo) and satisfy v{0) = and v'(t) ^ vq > 
for all t ^ 0. As before, we denote by Wt a standard Brownian motion. Consider the 
Gaussian process 



Zt = CT [ ^/Ti^dWs, (A.l) 
^0 



'0 

whose variance is a'^v{t). We will consider zt as a Markov process and introduce the 
notation F^''-^{zt G •} = P{zt ^ '\zs = x}, t > s, for its transition probabilities. Their 
densities are given by 

y^fit,y\s,x) := ^ y} 

= ^ ^ ^-iy-^)V^'rMt,s) (A 2) 

o- y/27rv{t, s) 

where v{t, s) = v{t) — v{s). 

Let d{t) be continuously differentiable on [0, oo) and satisfy d(0) > 0. The object 
of the "level-crossing problem" is to determine the density i/j{t) of the first-passage time 
T = inf{s > 0: Zg > d{s)} (which wc will call "first-passage density of Zf to (i(t)"). This 
problem has for instance been studied in [Dur85, RSS84, Dur92, ST96]. The aim of this 
appendix is to establish expressions for ip^t) useful in our particular situation. 

As D. Williams in the appendix to [Dur92], we will use two integral equations satisfied 
by ip{t). Let F{t) = f(t,d{t)\0,0) denote the value of the density of Zf at d(t) and let 
F{t\s) = f{t,d{t)\s,d{s)) denote the transition density at y = d{t) for paths starting at 
time s in a; = d{s). The Markov property enables us to write 

F{t) = f F{t\s)il;{s)ds. (A.3) 

^0 
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The second integral equation satisfied by ijj{t) is 



.o(i)F(t)- ['b{t,s)F{t\s)^l;{s)ds, 
Jo 



where 



boit) = v\t) 
b{t,s) =v'{t) 



d{t) _ d/ify 
v(t) ~ v'{t) 
d{t,s) _ d'{t) 
y{t,s) ~ v'{t) 



, , d (d(t) 



d fd{t,s] 



dt \ v{t, s] 



(A.4) 

(A.5) 
(A.6) 



with d{t,s) = d{t) — d{s). In the particular case of a standard Brownian motion, that 
is for (7 = 1 and v{t) = t, Equation (A.4) has been established in [Dur92, Appendix 
by D. Williams]. The general case is easily obtained from the fact that zt = crW^,(j) in 
distribution. 

Equation (A.4) suggests that the first-passage density can be written in the form 



1 



tp{t) = -c{t)e 



-d{t)^/2a^v{t) 



a 



(A.7) 



where c(t) is a sub exponential prefactor. In fact, the following bound on c(t) follows 
immediately from (A. 3) and (A.4). 



Lemma A.l. Let co{t) = bo{t) / ^/2^^v(t) . Then (A.7) holds with 

1 



it)-co{t)\ ^ 



Remark A.2. Note that Lemma A.l does not require v'{t) to be bounded away from 
zero as t varies. 

If, for instance, v{t) and d{t) are twice continuously differentiable, then s i— b{t, s) is 
easily seen to be bounded, and thus c{t) behaves like v{t)~^^'^ near t = 0. 

In [Dur92], an expansion of c(t) is constructed, which is shown to converge for all 
times t, under a convexity assumption on the boundary d(t). Taking advantage of the 
fact that cr is a small parameter, we can control the convergence of this expansion under 
a milder assumption on d{t), on a finite, but exponentially long time interval. Writing 
c{t,s) = b{t, s)/ \/2'Kv{t, s), we see from (A.4) that c{t) must be a fixed point of the 
operator 



{Tc){t) = co{t) - - f c{t, s)c(s) e-''^'''^/^"" ds, 
O" Jo 



where 



r{t,s) 



dis) 



d{tf ^d{t,sf 



v{t)v{s) 



v{s) v{t) v{t,s) v{t,s) 



d{s) d{t) 
v[s) v{t) 



1 2 



(A.9) 



(A.10) 



Remark A. 3. The exponent r{t, s) is nonnegative and vanishes for s = t. If r(t, s) does 
not vanish anywhere else, then the main contribution to the integral in (A.9) comes from 
s close to t. In the generic case dsr{t,t) ^ 0, the integral is at most of order cr^. If the 
functions involved are sufficiently smooth, one easily sees that the integral is of order a^. 
If r{t,s) vanishes in a quadratic minimum in s = i or elsewhere, then the integral is at 
most of order a. 
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It is well known that the most probable path reaching z at time t is represented by 
a straight line in the (u, 2;)-plane. Thus r(i, s) vanishes for some s ^ t \i and only if the 
most probable path reaching d(i) as already reached d{s). In that case, there exists a time 
u G (0, t) such that the tangent to the curve (f (s), d(s))s^o ^it (y{u)^d{u)) goes through the 
origin, i. e., {d{u)/v{u)y = 0. This situation can be excluded under a convexity assumption 
on d{t), which is equivalent to Hypothesis H3. 

The following lemma establishes the existence and some properties of a fixed point 
of (A. 9) under hypotheses tailored to our situation. We will employ it in Section 4.1 with 
v{t) = Jo* e^^W g{sf ds and d{t) = (2 - Si) e'^^'\ 

Lemma A.4. Assume that there are constants A, Mi,M2 > such that the conditions 

dis)v'is) - v{s)d'{s) ^ Av'{s){l + v^), (A.ll) 
\c{t,s)\^Mi, (A.12) 
M2v'{s) ^ 1 + v{s), (A.13) 



hold for all ^ s ^ t. Then (A. 7) holds with a pref actor c{t) satisfying 



W) - co(t)| ^ \' sup 

1 - £ vityi^ o^s^t 



whenever 



£ := 2Mi 



,-A2/4(r2 



■t + 



l + f(s)3/2 



co(s) 



A2 ) 



< 1. 



(A.14) 



(A.15) 



Proof: We shall prove that T is a contraction on the Banach space X of continuous 
functions c : [0,t] — > [0, oo), equipped with the norm 



|c|| = sup 



;(s)3/2 



As) 



1 + i;(s)3/2 

For any two functions ci, C2 G X, we have by (A.12) 



\Tc2{t)-Tci{t)\ ^ ||C2 -cil 



v{sfl^ 



(A.16) 



(A.17) 



Using Assumption (A.ll), we obtain 



d{s)__d(^_ r 

v{s) v{t) Js 



^ A 



d{u)v'{u) — v{u)d'{u) 

v{t,s) ^ ^ VW) - 



+ 2 



and thus 

r{t,s) ^ A^ 



v{t)v{s) 



+ 4 



_v{t)v{s) 

i^/W)-^Ms)? 
v{t, s) 



^ A 



2v{t,s) 



1 + 



v{t) I v{s) 



(A.18) 



(A.19) 



For the sake of brevity, we restrict our attention to the case v{t) > 2. We split the integral 
in (A.17) at times si and S2 defined by v{si) = 1 and v{s2) = v{t)/2. By (A.13), the 
integral on the first interval is bounded by 



)3/2 



4M2C7 



A2/40-2 Vy A2 



(A.20) 
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The second part of the integral is smaher than 2te~ 

-A2/4a2 because r{t, s) ^ for 

si < s < S2, while the last part is bounded by 



J S2 '^IV J S2 ^ 



This shows that T is a contraction with contraction constant e, and the result follows by 
bounding ||c — co|| = ||7""c — 7"0|| by a geometric series. Here denotes the function which 
is zero everywhere. □ 

Corollary A. 5. Let the assumptions of Lemma A. 4 be satisfied and assume in addition 
that there exists a constant M3 > such that 

d{t)v'{t) - v{t)d'{t) ^ Mail + v{tf/^) for all t ^ 0. (A.22) 

Then 



co{t) 



e M2M3 



^ c{t) ^ co{t) 



£ M2MS 

1 + ^ 

1 - £ A 



(A.23) 



1 -e A 

holds for all t > such that e = e{t) < 1, where e is defined by (A. 15). 
Proof: The proof follows directly from the bounds (A.22) and 

l + v{tl^J_ 
v{t)y^ coit) ^ A ' 

the latter being a consequence of (A. 11) and (A.13). □ 
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